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Abstract 

It has been shown that calculations based on a hadronic rescattering model agree reasonably 
well with experimental results from RHIC Au+Au collisions. Because of the large particle densities 
intrinsically present at the early time steps of Monte Carlo calculations attempting to model 
RHIC collisions undesirable artifacts resulting in non-causality may be present. The effects of such 
artifacts on observables calculated from the rescattering model are studied in the present work in 
two ways: 1) varying the time step and 2) using the subdivision method. It is shown that although 
non-causal artifacts are present in the rescattering model they have no appreciable effects on the 
calculated observables, thus strengthing the confidence in the results of this rescattering model for 
RHIC energies. 
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I. INTRODUCTION 



Calculations from a simple hadronic rescattering model have been shown to agree rea- 
sonably well with experimental results obtained from the Relativistic Heavy Ion Collider 
(RHIC) with v/i = 130 GeV and 200 GeV Au + Au collisions [J Q, [J More specifically 
these calculations provide good representations of the data for 1) the particle mass depen- 
dence of the my distribution slope parameters (i.e. radial flow), 2) the pr, particle mass, 
and pseudorapidity dependences of the elliptic flow, and 3) the pt, centrality, and azimuthal 
angle dependences of two-pion Hanbury-Brown-Twiss (HBT) measurements. Although the 
agreement between the model and data is sometimes more qualitative than quantitative, this 
agreement with such a broad range of RHIC experimental observables is still a noteworthy 
accomplishment for a single, simple, model. 

Many hadrons are initially produced in a relatively small volume in RHIC-type collisions, 
particularly at early times in the interaction. It is a challenging task for a Monte Carlo 
calculation of this type to deal accurately with binary collisions between particles in such 
a large particle density environment without introducing numerical artifacts which may 
affect the results. One such undesirable artifact that can occur when the interaction range 
between two particles is much greater than the scattering mean-free-path is non-causality of 
the collisions resulting in superluminal particle velocities jj, 0]. The method of subdivision 
can be used to minimize these artifacts associated with high particle densityjy]. 

The effects of such artifacts on observables calculated from the rescattering model are 
studied in the present work in two ways: 1) varying the time step and 2) using the subdi- 
vision method. The goal will be to determine whether non-causal artifacts are present in 
the rescattering model and, if so, whether they have appreciable effects on the calculated 
observables. This will provide a test of whether or not one can have confidence in the re- 
sults of this rescattering model for RHIC energies. Section II will describe the calculational 
methods used and Section III with give the results of the study. 
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II. CALCULATIONAL METHODS 



A. Hadronic rescattering calculation 

A brief description of the rescattering model calculational method is given below. The 
method used is similar to that used in previous calculations for lower CERN Super Proton 



Synchrotron (SPS) energies [7|. Rescattering is simulated with a semi-classical Monte Carlo 
calculation which assumes strong binary collisions between hadrons. The Monte Carlo cal- 
culation is carried out in three stages: 1) initialization and hadronization, 2) rescattering 
and freeze out, and 3) calculation of experimental observables. Relativistic kinematics is 
used throughout. All calculations are made to simulate RHIC-energy Au+Au collisions. 

The hadronization model employs simple parameterizations to describe the initial mo- 
menta and space-time of the hadrons similar to that used by Herrmann and Bertsch 0. 
The initial momenta are assumed to follow a thermal transverse (perpendicular to the beam 
direction) momentum distribution for all particles, 

1 dN m T 



ttit dm T exp (^) ± 1 

where ttit = \Jpt 2 + m o 2 is the transverse mass, is the transverse momentum, mo is 
the particle rest mass, C is a normalization constant, and T is the initial "temperature 
parameter" of the system, and a gaussian rapidity distribution for mesons, 

dN/dy = Dex P [-^^} (2) 

where y = 0.5 In [(E + p z )/(E — p z )\ is the rapidity, E is the particle energy, p z is the 
longitudinal (along the beam direction) momentum, D is a normalization constant, y$ is the 
central rapidity value (mid-rapidity), and a y is the rapidity width. Two rapidity distributions 
for baryons have been tried: 1) flat and then falling off near beam rapidity and 2) peaked at 
central rapidity and falling off until beam rapidity. Both baryon distributions give about the 
same results. The initial space-time of the hadrons for b = fm (i.e. zero impact parameter 
or central collisions) is parameterized as having cylindrical symmetry with respect to the 
beam axis. The transverse particle density dependence is assumed to be that of a projected 
uniform sphere of radius equal to the projectile radius, R (R = tqA 1 ^ 3 , where r = 1.12 
fm and A is the atomic mass number of the projectile). For b > (non-central collisions) 
the transverse particle density is that of overlapping projected spheres whose centers are 



separated by a distance b. The particle multiplicities for b > are scaled from the 6 = 
values by the ratio of the overlap volume to the volume of the projectile. The longitudinal 
particle hadronization position (zhad) and time (thad) are determined by the relativistic 
equations |9(, 

Zhad = r had sinh y (3) 

thad = T h ad COsh?/ 

where y is the particle rapidity and Thad is the hadronization proper time. Thus, apart 
from particle multiplicities, the hadronization model has three free parameters to extract 
from experiment: a y , T and Thad- The hadrons included in the calculation are pions, kaons, 
nucleons and lambdas (ir, K, N, and A), and the p, lu, rj, rf, <ft, A, and K* resonances. For 
simplicity, the calculation is isospin averaged (e.g. no distinction is made among a tt + , it , 
and 7r~). Resonances are present at hadronization and also can be produced as a result of 
rescattering. Initial resonance multiplicity fractions are taken from Herrmann and Bertsch 
sj], who extracted results from the HELIOS experiment ^(|. The initial resonance fractions 
used in the present calculations are: rj/n = 0.05, p/n = 0.1, p/uj = 3, 4>/(p + oj) = 0.12, 
rj'/r] = K*/u = 1 and, for simplicity, A/N = 0. 

The second stage in the calculation is rescattering which finishes with the freeze out 
and decay of all particles. Starting from the initial stage (t = fm/c), the positions of all 
particles are allowed to evolve in time in small time steps (normally At = 0.1 fm/c is used) 
according to their initial momenta. At each time step each particle is checked to see a) if it 
decays, and b) if it is sufficiently close to another particle to scatter with it. Isospin-averaged 
s-wave and p-wave cross sections for meson scattering are obtained from Prakash et al. [11.]. 
The calculation is carried out to 100 fm/c, although most of the rescattering finishes by 
about 30 fm/c. The rescattering calculation is described in more detail elsewhere []]]. 

Calculations are carried out assuming initial parameter values and particle multiplici- 
ties for each type of particle. In the last stage of the calculation, the freeze-out and decay 
momenta and space-times are used to produce observables such as pion, kaon, and nucleon 
multiplicities and transverse momentum and rapidity distributions. The values of the initial 
parameters of the calculation and multiplicities are constrained to give observables which 
agree with available measured hadronic observables. As a cross-check on this, the total ki- 
netic energy from the calculation is determined and compared with the RHIC center of mass 



energy to see that they are in reasonable agreement. Particle multiplicities were estimated 
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from the charged hadron multiplicity measurements of the RHIC PHOBOS experiment 
Calculations were carried out using isospin-summed events containing at freezeout for cen- 
tral collisions (6 = fm) about 5000 pions, 500 kaons, and 650 nucleons (A's were decayed). 
The hadronization model parameters used were T = 300 MeV, a y =2A, and 17^=1 fm/c. 
It is interesing to note that the same value of Thad was required in a previous rescattering 
calculation to successfully describe results from SPS Pb+Pb collisions []]]. 

As an indicator of the presence of superluminal artifacts, the transverse signal propagation 
velocity, v*t, is calculated for all particle pairs in the rescattering calculation for each time 



step from 



tpres tprev 



(4) 



VsT = \/V 2 sx + V 2 y 



where v s is the signal propagation velocity vector, ri(t pres ) is the position of particle 1 for 
the present collision with particle 2 at time t pres , r 2 (tp ret; ) is the position of particle 2 for 
its previous collision at time t prev , and v sx and v sy are the perpendicular components of v s 
in the transverse plane. Note that even though the rescattering calculations are carried out 
with relativistic kinematics, it is still possible for v s t to be greater than c for some particle 
pairs due to the details of how the two-particle collisions are implemented. 

B. Subdivision method 

The method of subdivision is based on the invariance of Monte Carlo particle-scattering 
calculations for a simultaneous decrease of the scattering cross sections by some factor, I, 
and increase of the particle density by /, where I is called the subdivision |4j. As I becomes 
sufficiently large, non-causal artifacts become insignificant. The present rescattering calcu- 
lation will be tested comparing pion observables from the "nominal" calculation, i.e. 1=1, 
with subdivisions of / = 5 and I = 8. Note that a subdivision of / = 5 has been shown to 



significantly suppress non-causal artifacts present in some Monte Carlo calculations 
Since the particle density increase is accomplished by increasing the particle number by a 
factor /, the computer CPU time taken per event increases by a factor I 2 . For the present 
study, 1640, 124 and 120 events were generated for the I = 1, I = 5, and I = 8 samples, 



respectively. Fortunately, the statistical value of the / > 1 events is I times greater per event 
than for 1 = 1. As an example, the CPU time taken to generate the I = 8 sample of events 
was 600 CPU-hours on 2.8 GHz PC processors at the Ohio Supercomputing Center. 



III. RESULTS 



Figures 1-8 present results from the hadronic rescattering model simulating hadrons pro- 
duced in RHIC-energy collisions for various time steps and subdivisions. Quantities plotted 
are v s t histograms including all particles in the calculation and several hadronic observables. 
Descriptions of how the observables are extracted from the rescattering calculation are given 
elsewhere {2! • All calculations are carried out for a mid-rapidity bin, —2 < y < 2, and for an 
impact parameter of 8 fm to simulate a medium non-central collision which should result 
in significant elliptic flow for the purposes of the present test. Statistical errors are shown 
either as error bars or are of the order of the marker size when error bars are not shown. 



A. Time step study 

Figures 1-4 show comparisons of various quantities for three different time steps in the 



rescattering calculations, At = 0.1,0.05 and 0.2 fm/c. The time step used in ca 



culating 
so the 



observables which have been compared with RHIC data is At = 0.1 fm/cQ, 0, [i 
other time steps used represent half of the usual step and twice the usual step. Figure 1 
shows v s t histograms. It is seen that v s t extends beyond c with an exponential tail for all 
time steps, indicating the presence of superluminal artifacts. For t> sT < c all time steps lie 
on top of each other whereas for v s t > c they differ somewhat according to what would be 
expected for the different time steps, i.e. smaller and thus more frequent time steps should 
enhance the superluminal artifacts, whereas larger and thus less frequent time steps should 
reduce the artifacts. For At = 0.05 fm/c an enhancement of ~ 10 — 20% in the artifacts 
over the nominal case is seen, while in the case of At = 0.2 fm/c a reduction in the artifacts 
by about 30% over the nominal time step is seen. 

The effects of these time steps on various hadronic observables are shown in Figures 2-4. 

Figure 2 shows ttit distributions for pions, kaons, and nucleons, Figure 3 shows the elliptic 
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FIG. 1: Transverse signal velocity distributions for all particles for various time steps. 



flow, V?, vs. pt for pions and Figure 4 shows two-pion HBT source parameters [3( vs. px- As 
seen in all three figures, even though the different time steps have noticeable effects on the 
artifacts (as seen in Figure 1), they do not significantly effect the hadronic observables, the 
observables for all three time steps agreeing within the statistical errors of the calculations. 



B. Subdivision study 

Since it could be argued that the time step study presented above may not reduce the 
superluminal artifacts enough to see significant effects on the hadronic observables, the 
subdivision test is also employed since it provides a greater ability to reduce these artifacts. 
This is seen in Figure 5 which shows v s t histograms for the three subdivision used: I = 1, 5 
and 8. While all three subdivisions mostly lie on top of each other for v s t < c, one sees 
factors of 5-10 reductions in the superluminal artifacts for I = 5 and 8. If these artifacts have 
a large influence on the values of the hadronic observables extracted from the rescattering 
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FIG. 2: Pion, kaon, and nucleon mr distributions for various time steps. 

model, it should surely be evident for the I = 5 and 8 cases. Figures 6-8 show the effects of 
these different subdivisions on pion distributions, elliptic flow and HBT, respectively. 

As was the case in the time step study, the hadronic observables shown in Figures 6- 
8 are not significantly effected by the different subdivisions used in spite of the fact that 
superluminal artifacts are greatly suppressed in the subdivision study. 

IV. DISCUSSION AND CONCLUSIONS 

It can be concluded from Figures f-8 that non-causal (superluminal) artifacts while being 
present in the hadronic rescattering code used in this study do not significantly effect the 
hadronic observables calculated from this code. It is possible to speculate on why this is 
the case. There are three main features of the code which might tend to reduce the effects 
of these artifacts: f) individual particles are allowed to scatter only once per time step, 
2) a "scattering time" of two nominal time steps is defined during which particles that 
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FIG. 3: Pion elliptic flow vs. px for various time steps. 



have scattered are not allowed to rescatter, and 3) once two particles have scattered with 
each other, they are not allowed to scatter with each other again in the calculation. The 
present study was performed for non-central collisions (6 = 8 fm), but it is expected that 
the same conclusion, i.e. that superluminal artifacts do not significantly effect the results of 
calculations with the rescattering code, would also be obtained for central collisions. This is 
because the particle density is not very different for central collisions in the model than for 
mid-peripheral collisions since the particle multiplicities are scaled by the overlap volume 
for b > 0, as mentioned earlier. 

In summary, two methods have been used to test the validity of the present rescattering 
model calculations in the presence of non-causal (superluminal) artifacts. It is found that 
the results of this model are not appreciably effected by such artifacts, thus strengthing the 
confidence in the results presented previously from this rescattering model for RHIC. 



9 



B 
E 

Q. 

C/) 

DC 



6 











7ut HBT from rescattering 
RHIC, b=8 fm, mid-y 



■ 



a 



0.2 



0.4 0.6 

P T (GeV/c) 




open: At=0.1 fm/c 
grey: At=0.05 fm/c 
black: At=0.2 fm/c 



♦ 







0.8 



FIG. 4: Two-pion HBT source parameters vs. px for various time steps. The ordinate scale for A 
is shown to the right. 
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FIG. 6: Pion pr distributions for I = 1, 5, and 8. 
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